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TEEOEY OF GROUND VTBRAJTIOl^S OF A TWO-BL.^J 3 E HSLICOPTBR 
ROTOR ON ANISOTROPIC FTEZTBLB SUPPORTS 
By Robert P. Coleman and Arnold M. Peingold 


SUMMARY 


An extension of previous worlc on the theory of self -excited 
mechanical oscillations of hinced rotor 118x108 has boon made. 

Previously published papers cover the cases of -tiiree or more rotor 
blades on elastic euppoirts (such as Inndins gear) having either equal 
or unequal support stiffness in different directions end liie case" of 
one or two blade rotors can supports having equal stiffness in all 
hox’izontal directions. Tho missing case of one or two blades can" 
unequal supports lias been ti'eated. 

Tho nKitheittatical ti’eai'anont of tills case is considerably loore 
con5)licated than ti.-.e other cases because of the occurrence of differ- 
ential equations with periodic coefficients. The character Is tic 
frequencies are obtained from an infinite -order determinant. 

Recurrence relatione and convergence factors are used in finding the 
roots of tho infinite determinant. 

Tao results show t 2 ie existence of ranges of rotational speed at 
which instability occurs (changed ecanewhat in position and extent) 
similar to those possessed by tho two-blade rotor on equal suppoi’ts . 

In addition;, the existence of an infinite number of Instability 
ranges which occurred at low rotor speeds and which did not occur 
in the cases previously treated is shovm. 

Simplifications occur in the analysis for the special cases of 
Infinite and zero stiffness in one of the axes. Tae case of infinite 
stiffness in one axis is also of special interest because it is 
mathematically equivalent to a counten-otating rotor system. A design 
chart for finding tiie position of the principal self-excited insta- 
bility range for the case of infinite support stiffness in one direction 
is included for the convenience of dosigaers. It is espocted that 
designers will be able to obtain sufficiently accurate information by 
considering orGy the cases of infinite and zero support stiffness along 
one directi CSX together with the cases treated previously. 
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It is kaown that rotatlng-wlng aircraft may escperiance violent 
vihratlons while the rotor is turning and the aircraft is on tno 
ground,. It has "been found tiiat tliese vihraticns can he oaplainod 
without considering aerodynamic effects and that ■fdaoy are due to 
mechanical coupling hotween horizcniial huh displacements and hlado 
oscillations In the plane of rotation. A theoi'eticeil analysis of 
this vlhra-blon problem Is given in references 1 and 2. Reference 1 
deals with rotors having three or more equal hladea on general 
supports and reference 2 deads with two-hlade rotors on supports 
having the samo stiffness in ail. dlrcctlonB . 

Althou^i in actual two-hlade rotary-wing aircraft, the stiffness 
of the supports along the longitudinal direction is certainly 
different from the lateral stiffness, tiie equality of the stiffnesses 
was assumed in I'eference 2 hecause it penn3.ttod the mathematical 
simplification of dealing with, differential equations having constant 
coofficients and it was helieved that a theory eirqjlcylng such an 
asstmoption would he sufficient to indicate the nature of the most 
violent types of ground inatahllity. 

Hie present paper gives a theoiretical investigation of the 
generaiL case of a two-hlade rotor mounted upon supports,, of unequal 
atiffu<iss along the two statianaiy principal axes. It thus 
generalizes the problem of reference 2, and rounds out the studios 
of ground resonance begun in reference 1. As was shown in rofo ranee 2 
a two-hlade rotor poeaesees different dynamic properties adeng and 
normal to the line of the hlados . Equations of motion with constant 
coofficients foi:* the problem hraated in reference 2 could he obtained 
by usiiy; a coordinate system rotating with the rotor. This procedure 
succeeded because the supports wore assumed isotropic (equal 
stlffnetss in all directions). 'When the supports are anisotropic, 
however, it is impossible to avoid the appearance of-porlodlc coef^ 
ficionts In the equations of motian. 

Ttc present method of solving the dlfforontial equations of 
motion follows closely the process oiigjloyod in i’oforonco 3 tor a 
v3.bratlon problem in two degrees of froedoa. The form of solution 
is oapiessed by an exponential factor times a courolex Fourier serlos. 
Substlt.ution of tho formal solution into the equations of motion 
yields an infinite sot of-algobraic equations and an infinl to-order 
determinant for the detorminaticn of the Fourlor coefficients and tho 
charactorlstlc frequencies. Hie subsequent analysis is concerned with 
methods of finding the roots of tho infinite determinant. 
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Becaiise the methods heroin er^^loyed may also be userul in 
other i’otor problems, particularly in those involving forward -fll^t 
effects, the mattiematical analysis is presented in same detail. 

It is eaqpoctod tliat dosigners will bo able to obtain sufficiently 
accurate ini'ormation by considering only the cases of infinite av zero 
support stiffness along one dirocticn together with the cases of 
references 1 and 2. In order to avoid the necessity for oxtensTvo" 
calculations, a; design-chart is included giving the location of the 
principail self -excited- instability range for iiie case of infinite . 
support stiffness in one direction. 


lERITATION OF TEE B®::\n}IONS OP MOTION 


The symbols lAsed herein are defined in appemdix A. 

The equations of motion are obtained from Lagrange *s equatiens 
and fi’om. tiie expressions for kinetic and potential energy. Four 
degrees of freedom of the x’otor system are considered: conpanents of 
deflection of the I'otor hub in the plane of rotation, "and hinge 
deflections of the two rotor blades about their vertical hinges. All 
motions are thus assumed to occtir in tiie plane of the rotor. TJie 
I’otor is assumed to rotate at a constant angular velocity a>. The 
analysis can bo applied to rotors without hinges by assuming an 
effective spring stiffness and hinge position to represent the elastic 
deflection of the blado. ■ . ■ 

The pertinent physical parameters are: 

a reicLial. position of vertical hinge 

b distance from vorticaJL Mnge to center of mass of blade 

Hjj mass of rotor blade . 

m effective mass of rotor supports 

r radius of gyration of blade about center of mass 

K^, spring constants of the rotor supports along tho 
X- and Y-directions , respectively 

Ep spring constant of hinge self-centering spring 
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I«t tlie origin of tiie X,Y- coordinate system "be placed at the 
iindisturbed position of the rotor hub. At time t eqxial to 0, the 
line through the blade hinges and rotor hub is assumed parallel to 
the X-_aals. After a time interval t, letHihe rotor hub deflection 
be z and hinge deflect! cais be and Pg# respectively, 'where z 

is the coicplex position coordinate laeas'ured in a coordinate system 
rotating with the rotor. (See fig. 1.) Tlien the poaltlonB of. tho 
centers of mass of the two blades, as measured in fiaed coordinates, 
will be , respectively, 

= Cz + a + bo^l^x) 

■ . - ‘ ■ (1) 

Z2 = (.2 - a - bo^^ 2 / e^^ 


Tie Idnetic energy of the rotor system 7 can be written as 

1 r r ^ 

,.T <zizi + zgzg + r^ Qa + >• 

i" . • 

+ g-m(z + i<'cz)(i - i<Dz) (2) 


The first term in equation (2) represents the kinetic energy of the 
rotor blades, including tho energy due to ro-fcation, and the second 
term ii3 the contribution of the rotor hub. 

U]pon e:xpandlng equation (l) into power series in and P 2 

( only small deflections from equilibrium being considered) and 
substituting into equation (2), there is obtained 

T w-gM{z + imz)(z - toz) 

-^Pl - P 2 ^(i^iA)z - io^z + £obz +, 

*01 - P 2 )C^iz + CDbz + cobz - biz} 

+ (b2 + r2)(p^2 ^ p^2^ _ ^^^2 ^ 322j~| (3) 

\ 

where only tlie terms ihat are quadratic in the variables have been 
retained, end M represents the total mass of the rotor sys-tem. 
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The potential energy of the system T 


is given "by 



. ^ (j-2e2iCDt 




w 


where K is the average stiffness end <iiK is a measure of the 
difference of the two principal Btlffnessesj that is. 


K F 

K =~ y _±^ 





As in refei-ences 1 and 2, siisplifications in the analysis are 
introduced hy replacing the hinge variahles Pi and P 2 hy the 
new variahles 9q and 6j_ 

% "Ksi + 3s) 

®1 = 1(^1 - %) 


In teiTBS of these new variahles -the expressions ( 3 ) and ik) 
heccme, respectively. 


T = 7^ ( z + icuz)(t - licz) + 


“b 




(5) 
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and 

TT ^ ^ 


(So? ♦ 9j_®) - ^ 52„-2tofr) (6) 

■By use of the Lagrangian form of the equations of motion 


dt 




the following equations of motion for the rotor system are 
finally obtained: 


(D + itu)2z + i|i(D H- ici3)%3^ 2 - = 0 

(D - jto.)2z - l(i(D - iaj)%3_ fp-Z ' ^ ze^^'*^ = 0 
(D + lxo)®z - (D - ia>)^E + 21 (1 + A^oj^ ° 


( 7 ) 

( 8 ) 

( 9 ) 


where t>ie notation 


D = 


dt 


has "been used; andL-the following couibinations of the original 
paramete>r3 hare been introduced; 



an , i-'« 

- 
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Equaticn (t) can be aolved independently of the o-toore since 
it is an equation in &q alone. The equivalent equation appeared 
also in references 1 and 2, and its solution represents in-phase motion 
of the "blades vitli no resultant reaction (except torsion) at the rotor 
hu"b. This motion will not he further considered herein. 

The problem is thus resolved into the solution of the three 
simultaneous equations (8) to (ID). It will bo noted that the terms 
with periodic coefficients in equations (8) and (9) disappear 

if = 0, that is, if Ky = S^. Equations (8) and (9) ar© tiius 
reduced to the problem treated in reference (2) . 


FOBM OF SOLIKION OF EQUATEOriS OF MOTION 


The eqxiatians of motion (eq-oatioias (3) to (lO)) are similar in 
mathematical properties to Mathieu*© equaticn, which occurs in tiie 
analysis of vibrating systesis of cno degree of freedom wito variable 
elasticity. (See reference k.) A generalized form of Mathieu’s 
equation was solved analytically by Eiil. (See reference 
pp. 413-417.) An extension of Eill's method has been applied in 
reference 3 to a problem involving two degrees of freedcm, snSTls. 
further development of the method of reference 3 ie followed in the 
present paper. 


Equations (8) to (lO) constitute a system of linear differential 
equations with periodic coefficients. Three second-o rde r .equations 
possess six linearly independent solutions that, according to the 
Ploquet theory (reference 5, P* 4-12), are of the form of an 
exponential fantor times a periodic function of tine. Particular 
solutions are of the form 

z = e^atp(t) + e-^a^t) 


2 = e^“^'*^Q(t) + o“^“^-^''T(t) (n) 

01 = e^^Vt) + e-^'^(t) 


whei ’0 is known as the characteristic exponent, and P(t), Q(t), 

and P(t) ere periodic functions of period 
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Since P(t), Q{t), and E(t) are periodic functions of t, 

they can he represented hy coij^ilox Fourier series, and equaticns (ll) 
bocama 




CO 


CO 




— CO 
CO 


(22ci>KDa)it X - -(2Zci>K%)it 


i = 'V Bj ^(22<o+Xia)it 


Z. 


^CO 
CO 


1 “ y_ 


(22oKUa)lt 


— 00 

- ,gO- 


-pCO 

JSSL. 


^ g-(2Ic>kUa)lt 

- -(22a>K2ia)it 

Cl e 




CO 


(.12) 


where Aj , Bj, end Cj are cc® 5 »lex constants. 


ICquations (ll) and. (12) show that, when the rotor system, is 
stable and oig^ is real, t}io motion not only is not simple harmonic 
as waej the case in refex'encos 1 and 2, but, in general, is not even 
periodic. The motion can bo said to consist of a fundamental 
frequency plus *'haimonics” of fi'equency cog^ + 2^05 wiiero I is 

any integer. From, equations (l2) it is seen that the value of 0 % 
is not uniquely dotermlnato, since <Sg + 2loi also satisfies 
equations (12). (The imaginary part of c% is definite, however.) 

It can be shown furthermore, that, corresponding to each value 
of 03;:^, -oua is also a soliition. Only those three values of c% 
for w>ilch the real parts lie between 0 and co need therefore be 
considered. These values will be referred to hereinafter as the throe 
**princ!ipal” values of These values of co^ differing from the 

principal value by 2l<a, or having opposite sign, will be referred 
to as "harmonics” of the corresponding prlncipial value. 


Eiince z has boon defined as a position coordinate in a 
rotating frame of reference, the values of o:^ can be interpreted 
as th^ natural frequencies of the rotor system in rotating coordinates. 


SOLUTION OF SITUATIONS OF I-K)TI0W 
Detoiminantal Equation 

If the formal solution (equation (12)) is ccmbinod with the 
equations of motion (equations (8) to (lO)), and the coofficionte of 
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each esponential tins factor is separately equated tb zero^ an 
infinite set of hcmogeneous equations is obtained. lEiese equatibrxs 
can "be separated into two independent sots. Each, equation of one 
set is the conjugate of an equation of iiie other set, and only' one 
set need ho considered. Thxxs 

Aj --^Sz+l 

- + (21 + = 0 (13) 


+ (22 - Bj - . 

+ (Xij^ + (21 - 1)^^ Cj = 0 ’ • (14) 


• jo)a,+(2Z + l)o^ + Ijng^ + (21 - l)o^^2 


+ 2ia + -(cJOa, + 2to) +A^o> +A2 = 0 (15) 


where I taices on all integral values frcaa -00 to 00. 

In order that the values of A^, Bj, and Cj not equal 
zero, the determinant of lixe coefficients of Aj, , Bj, and Cj 
iruiBt he zero. This determinantal eguaticai is - . 


A 
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a.4,4 a.4^_3 

0 0 

0 

0 

0 

0 

0 . 

^-3,-4 ^-3,-3 

^- 3,-2 ° 

0 

0 

0 

0 

0 . 

0_ a_2,-3 

^■-2,-2 ®'-2,-l 

0 

0 

0 

o' 

0 . 

0 0 


*-1,0 

0 

0 

0 

0 . 

o'“" 0 

^ ®0,-l 

®tD,0 

*0,1 

0 

0 

0 . 

0 0 

0 0 

*1,0 

*1,1 

*1,2 

0 

0 . 

0 0 

0 0 

0 

*2,1 

*2,2 

®2,3 

0 . 

o 

o 

0' 0 

0 

0 

^3,2 

^3,3 

‘^,4- 

0 0 

« • 

0 0 

• « 

• 

0 

• 

0 

» 

^4,3 

• 

*4,4- 

• 


= A(OJg^) =« 0 


(16) 


where 


A .4 = -X * 

W ” 3cu)2 

®0,0 

-4,-3 “ -^3 

®0,1 

_,/'c0a - 30)^2 



*1,0 

Ai<j 32 + Aa 


■3^-3 „ 2a>)^ 



A^o)^ + Ag 


V 


V 

2 


kM 


(c% + a>)2 


-■ b 0 
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a 

-3,-2 

^-2,-3 

^- 2 ,-?. 

®--2,-l 

®--l,-2 

®--l,-l 




= -1 + 


iz/u 

(o)a - 0i)2 


A 

(aig. - o)2 


A g/M 


(ojg, - cd)2 

.1 . 

K - 0)2 


^- 1,0 “A 3 



= A K/M 

1,2 (o:^ + <n)2 


A K/m 

^,1 (lll>g^ + tu)2 




"•4,4 


= -1 + 


K./M 


(i% + 2m)2 


The determinant has "been someviiat simplified l>y multiplying and 
dividing the rows and coluiims hy variaus quantities, and the 


pair ame ter 



has "been suhstituted for its equivalent 


U 

2(1 + r2/ij2) ■ 


Let tiiis infinite determinant he • T'l® prohlem consists 

in solving the equation A(o^) = 0 for its roots ca^. These roots 

will he Infinite in number, consisting, of the three prl ncj^pa l veilueB_ 
of co^ , plus al l the ir harmonics . The values of _ , as a 

function of cu/ \J~sJyi, are seen to depend only on tlie values of the 
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three nondlmansional parameters A]_> 
stiffjiess ratio parameter ^./K. 


A2 

k/m' 


and A^j ariA the 


.'i determinant of infinite order has meaning? caaly insofar as 
it is defined as the limit of a determinant of finite order. 

Define as the determinant of order 6n -3 formed from a 

squarf) array of A(a)« ) centered on the term -1 + . 

■ 

a 

This termi; which oi’l{?lnally was associated with in equation (ip)/ 

will ho referred to hereinafter as the "origin” of the infinite 
determinant A(cua,)* choice of tills term as center of A{c%) is 

purely arhitrfiry, and it was selected solely for reasans of symmotr’y. 


Thus 


A(i%) “ Ai(^) 

n->c« (17) 

The limiting values of the roots of the equation 
- «= 0 

as n hecomes Infinite will he tiio values of the roots of the 
equatlcaa 



■». 


A(o)a) => 0 


The inethod of calculating the roots of A(c%) = 0, hy 
successively calculating tiie roots of d^(coa) = 0 for larger 
values of n, is entirely too tedious. Instead^ the method of 
reference 3 will he followed. This method involves tlio calculation 
of the value of ACana) for several- specified values of cUg^. The 
roots of A(c%) = 0 can then he obtained from a trlgoncanetrlc 
equation Involving ttie roots .end the calculated values of A(o%). 


Auxiliary Doterminants and Recurrence Eolations 
for Calculating A(r 0 g) 

Before the trigonometric equation is derived; it is convenient 
to have a systematic ntunerical procedure for determining the value 
of A(c%). As n heocmes infinite; the terms of extend to 

infinity both above and below the origin. By expa n ding d^(£%) in 
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terms of the elements of the column containing the origin, it can 
he expressed in terms of anxi3JLary deteiminants that extend to 
infinity in only one direction. Recurrence relations can then he 
obtained that give the value of these auxiliary determinants. 

fhe auxiliary detennlnanta are minors of 
defliied as foilows: 


Cjj^(o3„ ) determinant of order 3n. -2 consisting of the terms 
below and to- iiie right of the origin; that is, 
determinant having first row and column beginning 

K/M 


with term -1 + 


+ cd)^ 



deteraiiiient 

omitting 


determinant 

omitting 


determinant 

omitting 


determinant 

omitting 


determinant 

omitting 


of order' 3h - 3 formed from C^(a^) 
last row and column 

of order 3 sl - ^ formed from 
laat~ row and column 

of order 3n- - 3 formed frcm Cp^(t%) 
first row end column 

of oxxLer 3h - ^ formed from I^(c%) 
last row and oolxnan 

of order 3n - 5 fomed from l^(cJDg^) 
last row and column 


ly 


by 



ly 


ty 


The following three determinants will also bo needed: 

Gj^(o!^) determinant of. order 3a - h formed from I^(c%) by 
omitting first row and column 

Hn(«%) determinant of order 3n - 5 formed fi'om 
omitting last row and column 


determinant of order 3a. - 6 formed from En(^a) ^7 
omitting last rCT? ond column 


Determinants similar to the foregoing can be formed in the 
same manner from the upper half of ^(%.) • Denote these determi- 
nants by the subscript -n Instead of n. It is seen, iiowover, that 
tlieii' values can be obtained frcm Hxe values of the determinants 
already defined from the lower half of by merely 

replacing with -cOq^ (for exan 5 >lo, C-nK) = <hi( -<%))• 



Eigianding ^(cua) In terms of tiie elements of tLie coluom ccntaining the origin gives 




/ 




- uNS / . Ak + 

(— 1 ,Iv(-iU^) f: fetv.^ I T > n \ 

^ y - nv-^x ' \^ OJ j 


; (IB) 

The aimiieiy determinants C^{o^), endiE^Ccpa) satisfy liie foUo^ 

recm-ronce relations (obtained by sxpan^Jng each in terms of the elements of its last rCTf)* 

"■ 1 r . : 12 

- A3 




-1 + 


^ + (2n ^ l)co]^ 


% + (2n - l).cD 


^ C %) = 


At^CoS +Ao 


% + (2n - 2)cD 
cua + ( 2 n - 3 )q) 


+ (2n - 2 ) 03 ]' 


|c% + (2n - 2 )cq 


SnW 


Cn-lW 


> 


VlC ^) 


- <-l + ^n~lW - njT 

1^ [o^a + (2n - 3)Gi]'"J t (2n - 3)t^ 

Tlia deteim^ts V^{^) , ^ and also the system ^ » 

satisfy the same recurrence relations as 0^(0:^), r^(ft^)j end 2^(00^);. 

3 ^Sp©ctiT©J^ • 


(19) 


H 


■'■ 1 


I I ‘ I - '• ■* 

I ' I : ■* 17 


It. ^ I > 
I I **! ''' - I * 

I 'If 


1J9TT ‘OH HI VOTH 
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The Tallies of these nine determinants can he found from tiae 
recui-renca relations (aquaticHiB (19)) and the following Initial 
Tallies, obtained directly from equation (16): 

K/M 


(zac/M )2 
CPa + 


= - 

% ^2 
(Pa. + 


r 


11 

(Pa + p)^-_ 




1 .. - 

(Pa + 

1 

K + 


o 

A2P‘~ ■‘‘As 

Coia + 



> 


(20) 


Go(ciig,) 


+Ao 

-1 +— ^ 

(cugi, + 2ui)^ 


r 

^ Aio)2 + A^ 


L - e/m 

(oig + a«;2 


iPa + 30i)^ 



By use of tJio initial conditians (equatlcns ( 20 )), the 
recurrence relations (equations (19)), and equation (iS) the 
Talue of A^(cSg^) can he calculated. The Talue of A(cUg^) will 
then he tlie limiting Talue of "becomes Infinite. 


The BehaTior of Ajj(cJ2.) ^or Large Values of n 

So far it has hoen tacitly assumed that the determinant A(a>g^) , 
as defined in equations (16) and (IT)^ is conTergent, and furtlior, 
that it remains a function of co^. in the limit as n hecomes 
infinite; that is, it is not Identically equal to zero, independent 
of the Talue of It will now ho shown that the function A^(o^) 

does hocoms zero in the limit, independent of hut that when ^ 

A^(o>g^) is dlTided hy an appropriate function of n, a now function 
Fj^(c%) will he obtained which will ho conTergent and remain an 
unambiguous function of in tiie limit. 
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^Ohe derivation of tbte appropriate function "by which to 
divide evidently depends upon tlie "behavior of ^(<^) 

as n "becomes vei'y large. As n becomes Infinite, the 
recurrence relations (equations (19)) beccane 


C^ = -D^-A^ ^ (eia) 

=• -1^ --^sCn-l (21b) 

Sn = -^-1 (21c) _ 

with identical equations for L^, and and for Hn, 

and Equations (Sl) ai*e readily solvable since they conetituto 

a systam of difference equations with constant coefficients. 'fi:ey 

are satisfied by solutions of the form 

Cn = OoiP- (S2a) 

where Cq is some arbitrary constant and Ic is a constant to be 
determ;.n0d. From equations (21c) and (21b), respectively, 


(22b) 

and ;i-. 

= Cq^I ( 22c) 

Combining equations (22a), (22b), and (22c) with equation (21a) 
and dividing througft by Cq}s?i -1 gives 

h = 2A3 - 1 (23) 

Thus, by use of equation (18), it is seen that for large values of n, 
'^^a) Since by definition A3 must have a val’^ 

between 0 and l/2, k2 must lie between 0 and 1. d3hue, in the li m .i 

lim 6n(Pa,) = 0 
n^ oa 

independent of co^ . 


Cf. - 
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Consider the function 


" 


j-2n 


The equation ~ 0 will ohviously have the some roots 

for cu^ as does = 0. 15ie function ^’nC^a) 

advantage, however, as is seen from, the preceding discussion, of 
remaining an unai±>iguous functicn of o>a_ in the limit as ti 
■beccmes infinite. Define this limit, as 


FQ^a) = lin Fn(5%) 

n-> CO 


= 11m 

n-v- “ 


k2n 


(2U) 


The primary problem can now he redefined an the problem of 
determining the roots, infinite in number and ccnslsting 
of <^a 2 ^ their harmonics , of tho equation 


F(oJa) » 0 


(25) 


Evaluation of Roots of Equatioa(2[>) 

Ihe following trlgoncsnetric expresBicn for E(ciJa) "wiH now 
be derived; 


■mf \ — -.j Ai(ta) 

E(c%) = lim 

n-^ 00 j5-2n 



(26) 
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The function F(aJa) is seen from equation (24) and 
equations (l6) to be periodic of period 2ui, to have roots 

± ^('^2 * *(^a 2 * vhere s is any 

integer, to have second— order poles at 0)^ = ±2soj, and to have 
fourth-order poles at u>g^ = ±(2s + l)oi. Liouville’s function 

theorem states that a function of a complex variable (in this 
case, 03^) that is analytic everywhere in the complex plane, 
inclU'iing the region ct infinity, must be a constant. It will be 
shown that F(<Ug) is finite at infinity (except if proceeds 

to infinity along the real axis). If the poles along the real 
axis could be eliminated by forming a suitable function of F(o)g^), 

without at the same time Introducing new poles, then that function, 
by Liouvllle's theorem, must be a constant, 

Such a function of F(cJa), which is analytic everywhere in 
the complex plane, is 



■where , c %2 ’ ‘“a 3 three principal values of o>a • 

The function thorofore a ccmstant. The value of J(cu^j 

found by making cOj, approach Infinity along the imaginary a::is is 


J(Po^ " 5’(oo) 


whero r(«») is the vaJ.ue of F(<Dq^) , as ■becomes Infinite . 

The value of F(co) cen bo found by letting ^ 

and then letting n-> w . Fi'’om 'tho form of f;^(co) it follows that 

f\^(co) = N2 jj(oo) = - 

The rocurronco rolations defining oquatioaas (2l). 

The expressions and may be ol.lnlnatod from oq-aatian.s (2l), 

whio}i ipivoB 

On = -{1 - %)Cn-i 
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Tlie initial canditions (equations (20)) reduce to C2_ = -1, from. 
» which it follows immediately that 

An(co) = - (l - 2.43)^'- 


=k2n-l 


from, equation (23) . 13 ierefore 


F(oo) = lim 

n-^ai 


Aa(c=.) 

t2n 


lim 
n — ^co 





(28) 


Equations (27) and ( 23 ) iamediately give equation (26). 

Aftei' equation (26) has "been obtained, the problem of 
determining ^2.^ ‘^'2^ considered theoretically 

ccnrplete, for equation ( 2 . 6 ) is really an identity in co^. Sxippose 
that cOg^ is assigned any specific valtie in equation (26) and 
is computed to a certain degree of accuracy. If these computations 
are made for two more values of ot^, all different, equation (26) 
will have yielded three equations in the tiuroe unkncvms ~ 

and ®J-9se equations can then be solved for the principal 

valuos of c%. Any degree of accxiracy nay be achieved by carrying 
out the computations for_. E(aijj) to a sufficiently large value of n. 

Ihe foi’e going procedure can bo systematized by rewriting 
equation (26) as - ^ 




-r-B 


. ■» - 4 -J 
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A ccn^'enlent choice for the three arhlbrary values of 
is tBgi = 0, (li, and (n/2. 

Kie o33)lioit definitions of K(o), K(l), and K(l/2.) tlien beccme 


K(0) . b 1„2 (^) BlnS 


e:(i) 


K(l/2) = 


1 - sin^ 




1 - sin2 


y r - -fel 


sin^l 


2 - 


2(H_21 


1 4 2 

— - sln‘== 


xao/ 


>(jo) 


The eauatione for evaluating K(o), K(l), and K(l/2) are 


K(0) = lim 


CO. 


sin-^ — ^ 


s(^ 

Van 




E(l) = iim 


cOg^— ^ CO 


kF(co^) 


cos^(^) 


> 


K(l/2) * lim 

0)/2 


iF(cOa.) 


(31) 


Carrying out the limiting operations indicated in equations (31)^ 
and using tho auxiliary detormlnahts 0 ^,( 1 %), and sc fori*, give 



\ 


fl 


% 


r 


K ( 0 ) = llm 

„ . O . . . Pn -1 

- ito“U - 2A . 3 ; 


> 




K ( l ) => llm 
n — > <» 


Jt 4 2 (< o ) 


Xb^l - 2 iiy 


^ ( 32 ) 


tr*/ m. "\ Jf-rn C 

IV \ i / r:y - UJ - Ui . 

n — ^00 


-1 +• 




. 8(1 - 2 i ^) 


2n-l 


y 


where tha quantities in or^ichots nro ccmLvaniuutljf rojrose-tsd hj 


. 


and. i ir(l/2)j^j respectively# The quantities K( are used in numerical ccaoiputatlons 
as approximations to the functions K( ) . 

The foraulaa (32) for K(0), K(l), and 'K(1/2) ocmverge slowly with increasing 

_ • 1 . s #1* rvP ‘l-VvA n"P 

n* 1510 corivorgence can oe spaeaea up grtJEiuLj 'j-*- 

cdivergence factoi’S used in reference 3- A convergonco factor for 1C( Is a function 

of n antiroachinw the limit 1 as n hecojnes infinite, which, when matlplled hy K( 


ro 
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gives _an S^lproSSioa which converges rapidly with increasing values 
of n to the values of K( ) . The details of the derivation of en 
EipP-T^oi'^ie-te canvex’gence factor for K(o)j^ will he found- in 
appendix B. Convergence factors for K{l)^ and z:(l/2)n ere 
derived, in a similar fashion, liha resultant expressions are 


where 


E(0) = Urn 

n — ^ rj^ 


K(0) sin^ ^ 
n 2 


J=1 \ 




E(1) * lim 

n— y oo| 


K(l/2) •■= 11m 
n — > « 


K(l) cos2 

n. g 

TT ^ ^ — ') 

\ (2.^ - 1)^/ 

K(i/2)j^ cos 
gn 

J=1 


fr(, . 

•I- \ (.S3 - \)V 


Q = 


T'* 

^ - A 3^ +A^cd^ +A 2 + 2h.y?- 

_ 2A3^ 


I 

> 


(33) 


F<sr a given value of n, the quantities in brackets eure found 
to bo better approximations to the rospoctivo values of K( ) than 
K( alone . 

The method of obtaining the values of CDg^ may be Bunaaarir.ed an 
follo?ra: by use of tiio Initial coiaditianB (equations (20)) end the 

recurrence relations (equations (l9))^ the values of the determinents 
^(^)^ Cr^(*<c)j Cj^('^/2), Cjj(-(o/2), I^jj(tu/2), end Lji( -^/?) 

cen be coE^iUted for increasing values of n. ¥lth the substitution 
of these values intor-equatlons (33)/ end with the use of equatlcns (32) 
a.pproxj,mate values of K(0), X(i), and K(l/e) can bo oanq^utod. 
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The process appears to he rapidly ccnvergent with n, especially 

for large values of (a/]fK/u. The values of end 

can then he found from equations (3l)> the definitlcsis of K(0)^ 
K(l), and K(l/2). 


Conditions for Stability 


From eqiaation ( 13 ) the cmditlon for stahllity of the system 
is seen to he that all three values of mst ho real numbers. 

If any one of them is coamplex or pui'e imaginary, then one of tlie 
terms in the solution (equations (ll)) will increase indefinitely 
with tiie time, tiie motlcxi therefore being unstable. . This ccaodition 


iirplies that the expressions sln^j 


and sin2 


\2taJ 



sin' 



all ai’s real positive nuniberB less than or equal 


to 1. The conditions for stahilitv can he eaprossod directly in 
teimis of K(0), E(i)^ and E(l/2; by moans of their definitions 

(equations ( 30 )). The throe equations ( 30 ) are formally equivalent 
to a single cubic equation 


x3 + h3^ + ox + d = 0 


/jtavA 

the roous X;j_, xg, and x^ of which are sin^^- -J and so forth, 

and the coefficients h, c, and d of -which are functicns 
of K(0), K(i), and E(l/2) where 

2 h = li-K(o) + 4K(1) - - 3 


= -6K(0) - 2K(1) + 8k(-|^) + 1 


2c 



2k 
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After some manipiilatiatis Inrolving tie. Descartes rule of signs, the 
necessary and sufficient conditions for stahility are found to he 


0 i 

-K(C) ^ 1 

0 1 

E(l)i 1 

Ji. CS 

8k(i/2) 1 


A = l3hcd - llh3d + h^c^ - l+c^') - 273.3 ^ 0 


The aumtity A is tho discriminant of the cuhlo equation. 


SPECIAL C/iSES OF GESSEAL TSEOEY 


Hireo special cases of the general theory are of interest. 

Tliese cases are the cases for vhich care, of tho principal stiffnesses 
Ky is reopeotlvelv' zero, equal, or Infinlle in magnitude in CQEjpariBcm 

with the second principal stiffness Ejj.. 

Case of Sy = Ex 

Tbe case of Ey = has heon treated in reference 2, If 
the equations (8) to (lO) reduco to the eouations of 
reference 2. In this speclel case, tlie motion of the rotor syatom 
hocomos simple harmonic., since all the coefficients Aj, Bj, 
and Cj In equations (12) ere identically zero except Aq, Dq, 
end Cq. 


Case of Ey = 0 

'J2ie special lind.tlng case of Ky = 0 is of interest in 

the case of a pylon of which the stiffness is negligible along 
one pr’inclpaJL direction withi interest centered on the frequencies 
involving the other principal stiffness. In the case of 
Ky = 0, tho function E(l) as given by equations ( 32 ) becciaes 
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identloally zero. lEiis roeult is also ©vldent fron the original 
definition of E(l) as given in equation ( 28 ), 'becanse one of 
tlie values of say is of necessity equal to +oi. (it 

■will "be reca3-lod tbat is the freqiiency as neasured in 

rotating coordinates. In fixed coordinates it wovild be zero.) 


It is possible to give inuch siir 5 )ler stability critericns for this 
case because there are only two E-functicns , K(0) and K(l/2), 

and two values of Uign ^nd <j 0 a 2 » determined. Ihe cev 

K-functiona may bo defined as follows: 



I'o « cos 



cos 



1 + ?js:(o) - 8 k(D 


(35) 


In terms of and the oriterione for stability boccsne 

0 ^ Kx 1 _(36a) 

0 ^K 2 1 1 (36b) 

-^+3/5^^! (36c) 

Given the values of Kx and £ 2 ? the values of anil 

c £) q ^ con be dotermined from equations (35)- A graph of the relati(m 

in equation ( 35 ) is given in figure g by means of which the real 
valt^es of end can be read off directly caace Kx and K 2 

are known, - 

A graph showing tJ.ie variation of and K, with <o/-yK^M, 
for the typical pereiaetero Ax = 0 . 1 , Ag^O, A-j^O.l, and 
K,.. = 0, is shown in figi\re 3 . By use of figure 2 the values of 
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and can te o'bts.ii'iod. These vcilues ere shewn in figure 4 

plotted against (a/'\J'SyjM.. CaloulatljcaiB are carried dewn 
, cu 

■CO = 0 . 5 . Ihe genoral^'behavior 'below this ppecd is 

discussed in the section entitled “CSeneral Behavior of Eotor 
SysteUi ae a Function of Eotor Speed.” 


Case of IC„ = CO 

The form.u3.as for the limiting case of Ky = 00 cannot bo 
obtained conveniently from the general theory. Instead of carrying 
out the limiting proconsj it appears pz’eforablo to boifln by treating 
the prob3.em as cno of only throe degrees of freedcan (two hinge - 
deflection coordinates and one hub-position coordinate i ), and 
by developing the theory along 3Lines similar to those used for the 
general ti'eataaient , In this way a system of two simultemooua equations 
witli periodic coofficienba is obtained, with the variables 6 ^ 

and X. Ihose oquatlone are solved in a mannor 6iml3.ai’— to that for 
tiie general case, the treatment being simpler, however, sinco tho 
solutl<m has only 'two principal values of 03|g^. 

Tlie dote.ils of tho solutlcai of tho equations of aDtion, to.gethcr 
with tlie final formulas for tho E-functicne, including ccnvorgcnoc 
factors, are given in appendix C. It is found that the some Kt 

and Kr,- occur os for Ky = 0. The oritarions for stability uxo 
exactly the same as those for Ky => 0 , tho ccsoditions of oquatione (36). 
Figure 2 can also bo used to detcimilno the values of from tho 
vali’.as of end Zo. 

A graph giving tlio veriaticin of , and with u>/Y Z^/hi for 

tho paraaacters A^_ = 0 . 1 , A 2 ^ 0, A .3 = 0 . 1 , and Ky. = w 1 b shown 

in figure p. In figure 6 the values of £%^/ 
ai’e shcim plotted against 


DISCUSSION OF KESULTS 
TypoB of Instability 


In|tablllty may occur as a result of the violation of any cue 
of tho igtability criterlc»i 3 of equations (i^). ‘Violatioia of o-r-ch * 

conditliDn is associated with a difforeait typo of instability, which 
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■would show up differently in -the niotlon of the ro-tor sys-oem. 
E:-:perienc3 with carapu-bs-tions indicates ^hows'rar , -tiiat -ere'cri tori one 
of most practical in^toi-bance for helicopters are 



-Z(0) ^ 0 
E(l) > 0 

Similarly j "the in^ortent criterions in tho lijaiiting cases of == 0 
and = CO ere 

. 

% >0 

Ko ^ 0 

If the condi'tion ^ ~ 0 Qdi* ■*" is wiola'bed, 

the o’thcr cenditions hein^ satisfied, then c%^ end ■will he 

conrolex conjusatos, and "blie rotor sye-ten will execute self-exci’bed 
■vihraticffis at fromiencios, in goneral, iricoEK&ensurete with ■the rotor 
speed. (Higher harmonics -will also he present.) This typo of 
instability ■will ho referred to hereinaf^ter as a "self -exci ■ted 
vibration." 

If the stability cr.ndltion -K(o) =0 (or K 3 _ = o) alone is 
violated, 'then caae of 'the ■values of ■will be a pure imaginary 

number. Physically, the rotor sys'tem ■vrill execute self -excited 
vibrations ha'ving a basic frequency, as seen in ro^tating coordina^tes , 
of zero. This behavior is similar to the ordinary critical -speed 
behavior of a shaft . Frequencies at higher harmonics 2ma w^ll 
also be present. This type of Ins^tability will be referred to as a 
" self -excited .whirling 

The third stability condition K(l) = 0 cannot be violated since 
K(i) as given by equation (32) cannot be negative. However, K(l) 
can be exactly equal to zero, (A similar s'tatement applies to K^-) 
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At such a point, "where the rotor system is cn "the hordor T"<r>i=t bo"bweon 
s-bahllLity and Ins-tahillty, csio of the veilues of e%, "Will ho ocuaJ. 
to itij. In fixed coordinates this result means that "the rotor sycton 

have a nattiral frequency equal to zero. Tlie rotor system "vd-ll, 
therefore, ho in I’eoonance "wl"th a steady force - a fox’ce of zex'o 
frequency. Hie aTnpIltude of "the zero-frequency term for- "the huh 
motion in such a si"buatipn can he shown to ho zero, hut "the hlodes 
"iiTill oscil3a.to. Also, hi^er harmonic terms, notably "the term of 
frequency Qx> (in ncnrotatlng coordine"tos ) , -will show up in t-io ht'h 
motion. Tb.le typo of vibration, which is a resonanco phenomenon and 
not a self ~exci"tcd "vibration, "will ha called a "s"teady- force rasonancG" 
vibration. 

lach of "the vibrations deacrihod - self-excl"tod vlhrationc, 
self -excited "whirling, and a steady-forco resonanco "vihratlctn - 
appeared in "the discussion of "the "two-hlado rotor cn equal supports 
(roforenco 2); however, thoro "the motlcna "Wore sirrolo harmonic, no 
higher harmonios being present. 


Gonoral Bohavloi’ of Rotor Sys"bcia as a 
Functlcn of Rotor Spoed 


T'le approxlma"bo location of "the ins"tahility regions Ciin 
easily ho found by examining "the limiting caco of * 0, that is, 

"tho case of zoro coupling ha"twoon "the hlado and huh motions. For 
sin^lioity, "the discussion is also restricted to "the caso of free 
hinges (A 2 = O) and Ky = «> . Tho Eq and K 2 functions hccamo 


% = 



cos 



\ 





E].imlnating "the rotor speed oj frem equations 


( 37 ) 


(37) givoe 
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Considei’ed as an equation in the vscrdables and Eg this 

equation represents a straight-ld.'ns segment (one of th.e lines in 
fig. 2) teininatod ty ’tihe K, and axes. Uie segjsant can "be 

shorn to be tangent to the curre =1. As (a decreases; 


the represcntati-ve point moves up end do'^n the lino segment; performing 
an infinite ntmber of ouch oscillatlans as co approaches zero. 

VJfcenever = 0, the point is at a solf-excitod-vhirling speed. 

“Ehe corresponding speed -is 




‘ 23 + 1 


where s represents any positive integer. Thus a seif— excited 
whirling will occur when the rotor speed is approximately equal to 1, 
l/3; 1/5, l/7; and so forth of the natural frequency of the 

hub V"^/M, Similarly it can be shown that there will be a steady- 

force resonance vibration whenever the rotor speed Is approximately 
equal to 1/2, l/4; I/6; and so forth of the hub natural 

frequency \jK^/hL, Finally self-excited vibrations will occur at 
rotor speeds approximately equal to 




2 s + 1 



Figure 7 shows the general pattern of rvjspcase frequency 
plotted against rotor spood for a small value of ii'.e mess-ratio 
paraiueter A-^. Tb.0 vai’iable 03/(0 rather then has 

been plotted ea oi’dlnate to avoid crowding of -tdie linos. Alcng the 
horizcn'fcal parts of the curves , blade motion prodcaninates over pylon 
motion. Pylon motion prodonrinates alcng tho slanting parts. 

Although the foregoing discussion was dovelcpod for the case 
of Ey = 00 , it is believed to apply equally well to the case 

of Ey = 0 and also to tho general case of E^ / E3- if the ro texr. 
hub is considered to hf-ve two natural froquoncios and 

oi'-ch frequency h:.ivixig associated with It an infinite set a£ iixfitability 
r.':inges located at approximately tho speeds given... - -- 
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Ccmparlsan. of Results for Different 
Values of 

Figures 4 and 6 give the princi pal v alues of 
plotted against the rotor speed (a/^Kx/n for Ky = 0 
and Ky = CO > respectively, both calculatecL -for the same set 
of paj’aineters At_, A2f and A 2* She calculations have been 

carried dcnm to 0 . 4 . The similarity between the two curves 

is strikl.ng. Sq far as tlae calculations have boon carried, each 
system shcrws tiie presence of one solf^xcited-vibraticn instability 
range, -one self -excited -wliirllng instability range, and. one eteedy- 
forco resonance speed A. If tne calculatlcns were carried to lowor 
values of o, further instability i-anges and stoady-forco resonance 
speeds, wotild appear . 


For coitrpariscn, the response frequencies of a two-blrAe rotor 
on equal supports (K^ = I^) for the same set of paraaetera is 

shown in figure 8. The fre<moncleB wex-e calculated frtaa the forEiula 
in reference 2 . Down to = 0 . 5 , this chai^t is very similar to 


figures 4 and 6. In ad-dition it ehowa one i-*ange of rotational speed 
a.t T/hich self -excited -vibration instability occurs, one range of 
rotational speed at which self -excited-whir ling instability occurs, 
and oni§" range of x^otaticnal speed at which a steady -force resonance 
speed occurs. Figure 8 differs px-incipally from the figores for 
Ky 7^ E;k; in that it shows no further instability ranges at low values 
of CO. • ' - - ■ ■ ' 


In referoncos 1 . .and 2 charts are presented giving the .location 
of the‘“ 8 elf-oxclted-vibraticn instability rcnge for varloxis values 
of the parameters An , A ^nd A-^* A similar chart for the case 

of a t'«ro-blade rotor, with K,, » oa la given in figure 9. In using 

*/ 

the chfirt, a. straight line is dx’awn ropi'esenting the variation 

p ''nCC^ +A„ 

of with the futxtlcn - The intex’sections of this 

... 

line with the appropriate A ^ curves give the beginning end end 

points of the Instability raaigo . The dashed line in figure 9 
illustx’atss tho mstliod for -the paremoters of figure 6 . 

Scpo obsex-vatlons concerning the relative location and extent 
of t!io various Instability ranges in figures 4, 6, and 8 appeal* to bo 
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applicable to a vide range of values of tie p?rajEstcrs h. 2 .) 

and A->. Thus the self -excited-vibration instability x'ango in the 
case of Zy = (fig. 8). is tfider (and hence tho vibration 
probably inore severe) tr.an tho corresponding ranges in the cases 
of ICy = 0 and = cd. (See figs. 4 and 6.) Also, tills insta- 
bility range occurs at lower rotor speeds in the case of ^ 
then it does in ti:e cases of =; 6xd = 0. The self- 

excited-whii'llng instability xx.nge is cansiderab3y" narrower 
for = 00 than, it is for the = S-jr case, and it is st3.ll 
narrovnx’ in the Ey = 0 case. 

In the general case of ^ location end extent of 

the inste-bllity rengC'S can be found falily ancuratoly by considoring 
the probleau as tho supei-pos’ltlon of twocrobltnas , ono of finding tlio 
Significant x'otor spcods rcfex-red to V G'B reference freouency 

with Ky assvjnGd infinite and the otlier of finding the sigrilficent 
rotor speeds referred to '^Ky/M es r*oforonco ffequ.onby Tfit:! 
assumed zero. With tiriQ foregoing discussion as a guide, s\ifficiently 
accurato design im'ona'tion can be obtedned withcr:.it extonsivo 
c’:.3.culationa for each val’io of Ey/K^ oncomterod in praetdeo. 


Effect of Hamping - 

■Mthou^ tlao of foot of dar^dng has not boon oxsninod mathe- 
mC'-ticoLly, because complications would, bo introduced in. the analysis, 
several inf oroncos fresa tlio darling investigations in roforencos 1 
and 2 can probably- ba safoly apiilied to tho rotor-system studios in 
tho pi’csent paper. Tlie nuxioroua instability r*angos occurring at low 
I’otor speeds, whiciT arc vci'y narrCT7 end reprosont a mild type of 
instability, arc probably cospletoly olinlnated by the prosonce of a 
Blight amount of damping in the rotoi' system. TIio prliuEiry solf- 
excl tod- vibration inst?.’.bility range can probably bo narrovred and 
oidminated by introducing sirt'ficient dainping into both tho rotor 
supports and tho blade hinges . 


APPLICATION TO imL BOTOPB 


It is easily shown tjiat the analysis for the case of . Ey = 00 
applies also to tiie case of a countorroteting I’otor system consisting 
of two equal two -blade rotors revolving at equal’ speeds and acting 
equally upon tlio same flexible mBiiber. TJxe rotbrs may be on the same 
shaft or on different shafts so Icng as the nonrotating flexible 
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member 1 b the same for hoth rotors. 
unoqti,al stiffness in iiio Z- and in 
the TJiideflectod bind© positions make 
stiffness axis. 


The supports, moreover, may h:>.ve 
the Y-directions, provided t^iat 
equal angles r/lth a principal 


The pi’oci. consists in showing that the energy expresslcsas for 
t^e dual-rotating system can ha separated into two independent sets 

which is of the same form as for a single rotor 
With E^ = cc. The resulting equationa of motion will thus also he 
the siame. 


:Eie separatiesn is acconqallshed hy introducing new vailahlos 

2 C^^oB " ®lneg) 

f t^os ■ ®°ueg) 

and 


Th == 


fOi_ 4 9 
V. -^OB 


^og) 


- — fe 

2 K 


+ 9 

'poe '^no 




where too subscripts pos and nog refei’ to the 6*3 defined for 

positive direction and for the rotca* tuiming 
in 1r.^^poslto direction, respectively. The energy ozprcsslons 

Dvicesns ^ JT- 


4 -^y^ + 


^ (' 33 ) 


■tr 2Kg f c. 
^ ” -b2 
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where M is the total mass of the Bystem. ' (M = m + 


Frcsa eguatlorLS (sS) it is seen, that £, is coupled only 
with. X and n is coupled on 3 .y with y. Equations (38) yield 
equations of 333oticai of iiie same form as equations (C2) and (C 3 ) . 


Tlie etahilitjr properties for the dual- rotating case are thus 
exactly the same as in the caso of Ey = 00 for the single-rotating 
two -h lade rotor. In parti cudai.-* , figure 9 can he used to find tiie 
location of the primary self- excited- vihratlon ins tahility” range. 


Tuo voluQ of A 
0,B A o = 


2 for the dual-rotating rotor is defined 

2mb _ mjj 

(m + 4iui3)(l + r^-ZhS) than o (m. + 2mb)(l ^ r^/h*^) 


as for the single -rotn ting rotor. All other parameters are the same 
for hoth cases. 


The quantity ^ sin oit can he interpreted physically as 

the x-compoD.ent of th.e displacoaent of tlie center of gravity of the 
h.lades due to hinge lef.lectioxis . quantity cos (at is tiue 

correspcaading y-oomyonent. 'Jlie aepai-atim of the variahles meansj 
physically, that tlrio 3 aotlon of the system can he separated into two 
independent modes, each of which involves linear motion of the 
supports along one of the prlncipeQ. stiffness axes. 

Siiidlarly, the stability of oounxerrotating rotor systems of 
six or more equal hlades can he detesmiined from the results of 
roferenco 1 with K.^ = 00. 


cowcniBioi^ 


The following conclusions are indicated hy the results of an 
iilvestigatlon of the prohiem of vihretion of a t\-o-hlade helicopter 
rotor on suppoz*ts that have different stiffnesses along the two 
principal stiffness •?o:es: 

1 . Many speed ranges are foiuid in which self -excited oscillations 
can occur- These oscillstions are of two types - self -excited 
vihration and self -excited whirliiig. There are a^o many speeds at 
which eteady-force-rescnance ■'.’dhi-ation may occur, 

2 . The self-excit-ed vihration, self -excited whirling, and steady - 
force resoni'nco speeds of higiiest rotor speed for each support natural 
frequency are recognized os corresponding to those of a two-hladc rotor 
on equal supports, hut chsinged somewhat in position and extent. 
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3. Mild, self -axel ted -whirling speed ranges exist at rotor speeds 

approximately l/3f l/5^ l/?* 6° forth of each support natural 

frequency. Steady-force resonance opeods exist at approximately 

1 / 2 f 1/6, and so forth of each support frequency. Self-excitod 

vibrations also occur at certain low rotor speeds. All these mild 
instability ranges are probably eliminated by the presence of 
medorate amounts of damping in the system. 

4. A familiai’lty with typical results of limiting cases of the . 

support- spring constants Ky =00, and Ky = 0 should 

enable a designer to avoid extensive calculations of cases of unequal 
support stiffness. In the general case of unequfil support stiffiioss, 
the location and extent of the Instability ranges can be found fadrly 
acouirately by considering the problem as the superposlticai oi' t^ro 
probloBis, one of fi nding significant rotor apeode referred to one 
suppoi’t frequency as reference frequency wi.th Ky assumod 

infinite and the other of finding the significant rotor speeds 
roferied to the other mipport frequency \i Ky/M' as reference 
frequency with asstmied zero. 

5 . Ihe analysis of a four-blade counterrotating rotor system in 
which the rotors cross along the principal etiffhees axes of the 
rotor supports leads to the seme oqmtions as thoeo considorod foi’ 
the special case of EL « <» , and the stability properties ore given, 
by the invostlgaticn of that special case. 


Lsngloy Memorial Aeronautical Laboratory 

Natlcaaal Advisory Cananittee for Aeronautics 
Langley Field, Va., July 22, 1946 
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APEEWLIX A 


SYMBOLS 


a 

®tDO’ ^l-» forth. 

3 ^, Cl 
h 

^<3- 

so fortii 
D 

F(o)g) = lim 


radial position of verticeQ. hinge 

elemsTita of detenainant defined in. equation. 

(16) 

Fourier coeffici emits in equation (l2) 

ccEplex conjugates of Aj, and Cj, 
respectively 

distance from vertical hinge to center of 
mass of "blsde 

minors of the d© terminate 
time -derivative operator 





A^CoJa) 


Jj n, s 

J(o>a) 

k = 5/I3 - 1 

Ki(D^/(X>) 

E(0), Ed), E(l/2) 
Ky 


integera 

fuiactlcn of F(<o„) defined by equation ( 27 ) 

Q< 

function of cOg/cc defined in equation ( 29 ) 

funotlona defined by the relations ( 30 ) 

functions defined by the relations (35) 

spring constants of the rotor supports along 
‘ the X- and Y-directicns , respectively 
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K 




L<K 




avereLge stiffneBS of rotor Bupporbs 


ITjj . spring constent of blade self- 

centering spring 

m — effective mass of rotor supports 

mass of rotor blade 

M total mass of two-blade rotor system. 

(m 4- £zi>ij) 

p(t),G,(t),R(t) periodic fvmctiona defined in equation (ll) 

P(t) ,^!(t) ,E(t) coEplex conjugates of P(t),Q(t), andR(t)> 

respectively 

Q . — . . f.Q|Tg-tant defined in e qua ti cans (33) 

r .... . .. .. radius of gyration of blade about center cxf 

mss 


t ...... . - time 

T kinetic energy 

■y potential energy 

X, y deflection of rotor hub neasured in 

X,Y-coordinate system 


X, Y fixed rectangular coordinate axes taJoon 

~ parallel to the principal stlfl’hoss 

direct! ons of the rotca’ hub 

z complex position coordinate of the rotor 

hub in rectanguJnr coordinate system 
rotating with angular velocity (zj, + lyp) 

s - ccmplex conjugate of z (Xj, - ly^,) 

Zf complex poslticai vector in X,Y- (nonrotating) 

coordinate system (x + iy) 
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l^l-’ ^2 


anga].ar hinge deflections of rotor 
hlades, respectiTely 


A 


A(0J^) 

«o " ?e) 

B, =|(Bi - Bs) 



^0 J 'HiJ ‘'To 


|i 


discrljQiii'?n.t of cubic eqtiaticfa 
:c3 + bx^ + cx + d = 0 

determinant of infinite order defined 
by equation (l6) 

determinant of order 6n - 3 formed 
from A(c%) 


blade Tariablps foi' counterrotating rotor 
mass ratio 



05 


constant angular Telocity of rotor 




cu, 


‘®l’ 


05 - 


characteristic esponent or natural, frequency 
of rotor system as Tiewed in coordinates 
/rotating at angiuLar Telocity co 

principal Talues of 03^ 
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APEENDIX B 

lEEITATrON OF lEE COFy/ERGENGE FACTOR FOE K(0) 

A convergence factor for K(0) is found "by finding a eln^iier 
function that changes with n In nearly the same way as K(0)j^. 

■Sien, if G denotes 11m G^, the expression 

n — >03 


K(0) a G 

^n . 

for a i^lven value of n is a "bet her approximation to K(0) than 
is K(0)j^ alone. A aui table form for G^^ is found from, a study 

of tJie be^iavior of K(o)j^ for inrge values of n. 

Tlie behavior of is studied by first observinc the 

behavior of Oj^(O) end l<n(0) for large values of n emd then 
inferring the behavior* of K(0)^ from equation (32). In the 

discussion of equations (22) it was shown that, as n becomes 
infinite, the ratio app3:*oach6S the value k (equation (22a)). 

A closed’ approximation to tlio value of this ratio can be written as 


w°) 

c„(0) 




(Bl) 


Whore p and Q 8d’e constants to be determined. 



K/.CA TN ITo. 11 84 


39 


Equatlcas (19) 13600(1116, for = 0, 




= 

W°) = 


-1 + 


-1 + 


-1 + 


E/t4 

(2n + 1)2 cd2 

h^2P^ + -A. 2 

(2n)^^_ 
K/M 

(2n - 1)^ 




^(0) - 


> (B2) 


cjo) 


where the second teim In the agnation for Ej^(£%) in equations (19) 

has been neglected as being of hi^.er order in p^^rs 

the terms retained. Eliminating E(0) and D(o) .ram the 

equations (B2) results in _ 


WQ) 

■ Cn(0)" 


1 IC/li 

-1 


— 

-1 . .. 

Q -(2n . 

(2n)^_ 


(2n - 1)'=’^ 



Upon expanding this expression into pwers of l/n and retaining 
terms up to and including those in l/n 


Cn.-!-l( 0 ) 

CjiC® ) 


-1 + 2^3 - i 


(2n)2 




. 2../^ - ) 


(B3) 


032 



IfO __ Np. llolt 

Ccffiipa^'ing equation (B3) with equaticn (B3-) 2?esults in 
P = 0 

^^1 -Ag) +A2 _q>^ + A 2 SA.gtt!^ 

^ C£)2(1 - 2A.3) 


Similarly 


I«(0) 




Therefore, from equations (32) 


E(0)n 



Sence, sn approximate value of 

g(0) - K(0)n-.l 

KCoj;- K(0)n ^(^)n+l 


K(0) 

,K(0) 


can he 
n±3- 


K ( 0)^^2 • • 


obtained fran 



The ri^it-hand side of equation (B^), which is seen to bo of the form 
G/Gj^, is a convergence factor. This ccnvergence factor ia the one 

used in equation (33) • 
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APESNDIZ C 


MATEEl-lATICAL ANALYSIS FOR THE CASS OF 




In teiiaa of the variahles 8 q , 9 t_^ 2^ ^ 'the espressiona 

fox' tho kinetic and potential enei'^ ane 


T = -^ + BHj 


Pic^Q-) sin <ot + cos cot^ 


^ -f ^ ^ I (^0^ ^ 


The three a(5uaticn.s of motion are 


^‘0 ^ ■^-'^2)^0 “ ° 

'• ^ d^ A . A 

X + — X - iJ. (0-, sin 'at) 


= 0 


~ ^2\ ^ cot + 9^ + +A^9^ = 0 

r " S') 


(Cl) 


(C2) 


(C 3 ) 


Equation (Cl) is identical wi -ax equation (T)- Equattcns (C 2 ) 
a.nd (03) constitute a STStoa. of tv?o linear BecorA-oi’dei' 
differential equations -witii periodic cocxfxclontc . 
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Equeiticns (G2) and (C3) are satisfied by solutions of the fona 



A^e 


i(c^+lm)t 


> 

if(%+(2+l)c3l^ 
■^2+1® “ 


(C4) 


where I talces on all odd Inte/^ral veiluss and and 

are constaats to he deterjnined along with the tJ^o principal 
Talueei of V®'! ^2^ ‘ constants Aj and 

in equations (C4) euro, of course, different from tiiose in (l2). 

Combining eqiTations (C4) with equations (C2) and (C3) and 
setting the coefficients of the various eaponontial time factccrs 
equal to zero 


1 + 


(o) + to)s 


u 


A, - 


21 


+ ^1+^ =0 


2i 




L" 


f trWy ■- j 



- 



- 2)co 

H~2 

“a + 

2 

(Da + (2 

- l)m 

a> -t- (2 - l)co 

H 


— 


•* — 

J 


> 


j +A 2 I 

. — - ■ 

[coa + - 1)“]^ I 


+ r<- 1 + .r 


2 Pz-1 ° 


Tho determinantal eqiiatian is then equal to 
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• • 

• "--3,-3 

^-3,-2 

0 0 

0 

0 

0 • • ■ 

• • 

• ®--2,-3 

a_2,_2 

^•-2,-1 ^ 

0 

0 

0 • V ■ 

9 • 

. 0 

®--l,-2 

®--i,-i ®--i,o 

0 

0 

0 • • • 

• • 

0 

0 

^0,-1 ®'0,0 

®-o,i 

0 

0 « • • 

• « 

0 

0 

Q ®-i,o 

®-l,l 

^■1,2 

0 • • • 

• ■ 

. 0 

0 

0 Q 

^•2,1 

®-2,2 

^•2,3 • * ‘ 

• • 

0 

0 

0 0 

0 

^3,2 

^3,3 * • * 


• 

• 

• • 

• 

• 

• 


9 

• 

« • 

# 

« 

* 


= A(ojg^) = 0 


where 

, Wm 

»->3 - c,^ . y^ye 




■op 


= -1 + 


■ ^1 g>'^ + 


^2 


_ 

^- 3 r 2 ” 2 




_ T +A2 

8'.2_2 - -- + (a^‘ . aij)2 


" ■*■(£%+ co )2 



-'^3 

ai 2 » 2 - 
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Kx/M 


“- 0,-1 ^ . co )2 


' ’ W6a m No'. iit4 

(oJa + 20i)^ 


As 

'^-10 2 


_ / a>a *• oA 

°-vi Kr^j 


The texTU 
detormJjient . 


(~' An jL.^O 

-1 + . Vill 


+ 3"’ 


(c% + ao)2 




Kx/M 


>= -1 + + 3u>)2 




I'c taken as the origin of the 


Define detonainan.t of order 4n-l formed 

hy taking a square array from centered cn tlie origin. Then 


A(a>a) = 11m 


n— >00 


Define auxiliary determinants from follovs: 






determ i nant of ordei' 2n.-l consisting of the terms to 
the ri^\t of and helow ttie origin term. 

determinant of order 2n-2 ohtained fream CJjj^(cDg^) 
hy omitting Its last row and column 

determinant of order 2n-2 obtained from 
by omitting its fii-st row and column 

determinant of order 2n-3 obtained from D^(uig^) 
by canittlng its first row and column 


T!a^> determinants Cj^(a)g^) and satisfy tiie follcwing 

recurronce relatians: 
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CnC^a) = ^ -1 + 


Kx/M 


[oj^ + (2n - l)<^^ 






P. 


cii„ + {Pn - l)tB 

8L 


ct> + (2n - 2>ai! 


^a-i. 


_iK) 


= 4 -1 + 


<o‘^Ai + Ao *~1 
<f _1 + ::a v 

i + (ai - 2)i^j^ I 

i_ -■ -!J 




(C5) 




As 


<D„ + (2n - 3)»|^ 

Og, + (£n - g )a>[ 




5die recurrence relations (equations (C5)) are also satisfied 
"by M^(cOg^) end •’ ^ replacing C end P, 

respectively. Tb.e ?nitial values are 


‘^l('%) ■= -1 + 
Do(<%,) = 

Np(03^) = -1 + 
Mo(sjp) = t 1 - 


Kj./M 


+ 

03)2 

r^M 

1 

+ 

0 ) 2 ] 

+Ao 

(ax^ + 

2o>)2 

+A 2 


au)2 


uj^A^ +A 2 


(cDa + ao)^ 


A 3 


1 - 


E^/M 


(cDa + 3a>)^ 


AoA 



•• 3^>y 


+ 2cd/ 
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E:qpajidlng Aj^(cDg^) in terms of the elements of the 
column containing the origin gives 


An(i%) = I -1 + 


A- 


+A 2 ' 




“a" 


r 


/“a - 


Aei n hecomes infinite , or as oig^ hocoaes infinite , the 
rocui'rence relations (equations (C5)) approach 

A 3 


Cn = Cn-1 


A 3 

■% = -Cn-i Vi 


(C6) 


Eqijiations (C6) are sabiefied hy a solution of the fom 


= ChP 










(CY) 


where h satisfies the eauation 




- - (i -A 2)4 “ 0 


(C3) 


The larger root of equation (C3) will he denoted hy k and tJie 
smaller root hy Although the complete solution of equs.t±an (G6) 

is of the form 

= Ckfi + C’k^^ 


for largo values of n the term in k^ hocomos negllglhle caaparod 
witli the term in k. 
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With the same values for k, and will have solutions 

similar to tliose of equations (CT)- ' Thus as n 1360011163 infinite 
Aq^(ci 3^) will vary as the quantity 3i^. 

Define the function 


F(o3a) 


Urn Fn(<ug,) 


= lim 
n— >03 


l£^ 


The function F(a>g,) is periodic in of period 2ca, has 

roots + 2303^ , t^osg^g ± 2 scd^ , and second-order poles at ^ 03 ^+ 003 ' 

for all integral values of s. Furfeermore j F(c3g^) approaches 
the limit 


E = lim F(o 3^) = lim Ejj^Coo) 

iOoL ^ ^ H— ^ 00 


*^2n(») 
« lim 

n — ^ 03 iphi 


k2.^ - k^ ' 


as cJ3^ hecaaes ix3finibe in a direction otlier then along the real ajds 


Foot the function 



(C9) 


The function J(tOg) is an analytic fmction of o>g^ everywhere . 
Hence, "by Liouville's theorem, ^ constant. By letting 

ca ^ — > + CD along the imaginary ajd.s, it is seen tiiat 1 ( 03 ^^ “ 
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S^'lDStitutlng -4E for J(^) Into equation (C9) results in 



Introducing 
the case K = 0 


E functicais defined siTidLleorly to those used for 
gives 



Carryi3ae; out the Hjait processes indicated in equations (CIO) filvos 
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^9 


= lim K_ 
n — >.co "Sa 


= Urn 
n — ^ 00 


+ A 2 )Cn^( 0 ) “ 3iAo^Cn( 0 )Nn( O ) 

-Jwj2 C2ii(<=>) 


ICg = lim 

n — i 00 ~ii 


= lim 

H- ^ CO 




-Vn 2 C^(oo) 


Finally, upon introducing appi'oprlate convergence factors, 
the Quantities needed in eq-jation (36) are 


Kj_ = lim 
n — > ■» 


K' 




sto 2 (f) 


. 1=1 ^ 




= lim <; 


coa-^ 


Pfra 


> 


n — ^ J3 




. 1=1 




( 2.5 - 1 )^ 


where 


p 2 = 


+ o^Vj_ + + co2(i + - Ic) 


Oi^Qc “ 112) 
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Figure 1.- Simplified meohanlcal system representing rotor. 
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Figure 4.« Frlnolpal ralues of plotted against rotor speed o) for oaee 
of A A 2 “ Of 3 “ 0*^t and Ky — 0, 
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Ti.gaT9 Oraph of and. Kg 

speed oi for = 0.1, Ag = 


as funotlons of the rotor 
O, Aj = 0.1, and K^. = oo. 
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Figure 6.~ PrlnoipAl Talaei of plotted against rotor speed (O for oaee 
of = 0.1, Ag = 0, A^ = 0.1, and = ®. 
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figure 6.- Responae frequenolea ot m two-blade rotor on ajmaetrlo eupporta for 
= 0.1, Ag =*0, Aj = 0.1, and Ky = 


■( < 
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